APPLICATION OF MULTIVARIATE SPLINES TO DISCRETE 

MATHEMATICS 



ZHIQIANG XU 



Abstract. Using methods developed in multivariate splines, we present an 
explicit formula for discrete truncated powers, which are defined as the number 
of non-negative integer solutions of linear Diophantine equations. We further 
use the formula to study some classical problems in discrete mathematics as 
follows. First, we extend the partition function of integers in number theory. 
Second, we exploit the relation between the relative volume of convex polytopes 
and multivariate truncated powers and give a simple proof for the volume 
formula for the Pitman-Stanley polytope. Third, an explicit formula for the 
Ehrhart quasi-polynomial is presented. 



1. Introduction 

Let M be an s x n integer matrix with columns mi, . . . , m n € Z s such that their 
convex hull does not contain the origin. For a given a € 1/ , consider the following 
system of linear Diophantine equations 

MP = a, (3e Z n . 

The number of non-negative integer solutions /3 for this system is denoted by t{pt\M ) 
and the resulting function t(-\M) on Z s is called a discrete truncated power. Discrete 
truncated powers, also called vector partition functions, have many applications in 
various mathematical areas including Algebraic Geometry, Representation Theory, 
Number Theory, Statistics and Randomized Algorithms. In general, one studies 
t(-\M) by generating functions and by means of algebraic geometry etc. [3J [3J [51 El 
128] . For example, when M is unimodular, that is when every nonsingular square 
submatrix of order s has determinant ±1, two algebraic algorithms for generating 
the explicit formula for t(-\M) are presented in [TH]. When s = 1, an explicit 
formula for t(-\M), which counts the integer solutions for the linear Diophantine 
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equation, is presented in [2]. Especially, Popoviciu [27] gave a beautiful formula for 
t(-\M), when M = (a, b) where a and b are relatively prime. 

In this paper, an explicit formula for t(-\M) is presented for any integer matrix 
M. In contrast to other ways, our method is based on multivariate spline functions 
and is inspired by the work of Dahmen, Micchelli and Jia [TU1 H21 EU [22], who 
exploited the relation between t(-\M) and multivariate splines, and demonstrated 
the piecewise structure of t(-\M). Moreover, we believe that the tool of multivariate 
splines, which have been developed in this latter theory, shed some light on problems 
concerning t(-\M). 

The main results in this paper are as follows. As the central result of the pa- 
per, an explicit formula for t{-\M) in terms of multivariate splines is presented. As 
applications of our formula, we first generalize the formula for partition functions 
in [2] and give a simple proof for Popoviciu's formula. Further, we show that the 
relative volume of a convex polytope agrees with the value of the multivariate trun- 
cated power at a certain point. We also present an efficient method for computing 
the volume of convex polytopes and re-prove the volume formula for the polytope 
related to empirical distributions, which is presented in [26] as a central result. 

The paper is organized as follows. To help make this paper self-contained, we 
introduce the multivariate truncated power and the box spline in Section 2. In 
Section 3, the discrete truncated power t(-\M) is introduced. An explicit formula 
for t(-\M) is described in Section 4. The proof of the formula is presented in Section 
5. In Section 6, we simplify the explicit formula for t(-\M) for some special matrices 
M. The explicit formula for partition functions in [5] follows directly from our 
formula. Section 7, containing two subsections, uses multivariate truncated powers 
to investigate the volume of convex polytopes and the Ehrhart quasi-polynomial. 
Particularly, in Subsection 7.1, an efficient method for computing the volume of 
convex polytopes is introduced, with re-proving the volume formula for the polytope 
related to empirical distributions which is the main result in [26 ; in Subsection 7.2, 
an explicit formula for the Ehrhart quasi-polynomial is presented. 
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It is necessary here to recall some previous notations (see [22] ). Throughout 
the paper, Z + and R+ denote the non-negative integer and non-negative real sets 
respectively. For given sets Dj, D 2 , let Ip^D^) = if D\ n D 2 — 0, otherwise 
let 1d x (D2) = 1. The linear space R s is equipped with the norm | • | given by 
\x\ = J2i<j< s \ x j\> wnere x — (xi,X2, ■ ■ ■ ,x s ) € M. s . Let A and B be two subsets 
of K s . Then A — B is the set of all elements in the form of a — b, where a e A 
and & € -B. The sets A + B and cA are defined analogously, where c € K. The set 
is the complement of B in A. A subset f2 of R s is called a cone if il + SI C fl 
and cJ7 C f2 for all c > 0. If a cone is also an open set, then we call it an open 
cone. Let Y be an s x n matrix. The linear span of Y, denoted by span(y), is the 
se t {J2 y eY a vV : fl y G M for all y}. The cone spanned by Y, denoted by cone(F), 
is the set {J2 y eY a yV '■ a y — ® f° r an V}- We denote by cone°(F) the relative 
interior of cone(Y). The convex hull of Y, denoted by [Y], is the set {X^gr a yV '■ 
a y > for all y and ^2 y£Y a y = !}■ Let M be an s x n matrix with integer 
columns mi, ... , m„. We denote by [[M]] the zonotope spanned by mi, ... , m„, 
i.e., [[M]] := "Ej=i a j m j : < a, < l,Vj}. Moreover, set [[M)) := -Ej=i a j m j : 
0<a, <l,Vj}. 

We shall use standard multiindex notations. Specifically, an element a € 
is called an s-index, and \a\ is called the length of a. Define z a := z" 1 • • • zf s 
for z = (zi,...,z s ) £ C s and a — (ax,...,a s ) G Z+. Also, set exp(ca) := 
(exp(cai), . . . , exp(ca s )), where c is a constant. For z = (zi,...,z s ) € C s and 
w = (wi, . . . , «j s ) € (C \ 0) s , set z/io := (zi/toi, . . . , z s /w s ). We denote by Pfc := 
Pfc(M s ) the linear space of polynomials in s real variables with total degree < k, 
where k £ Z. If fc is a negative integer, then we interpret as the trivial linear 
space {0}. Moreover, we also set P(R S ) := Ufco P fc( RS )- 
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2. Multivariate truncated powers and box splines 

Let M be an s x n real matrix with rank(M) = s. Throughout this section we 
always assume that the convex hull of M does not contain the origin. The multi- 
variate truncated power T(-\M) associated with M, introduced firstly by Dahmen 
8J, is the distribution given by the rule 

(2.1) / T(x\M)4>(x)dx = [ cj)(Mu)du, <j> e @(R S ), 

where £)(M. S ) is the space of test functions on R s , i.e., the space of all compactly 
supported and infinitely differentiable functions on R s . From (|2.ip , we see that the 
support of T(-\M) is cone(M). In fact, T(-\M) is identified with the function (see 

BSI pp-i2) 

2.2 T x M) = — " v V -, 

y/\det(MM T )\ 

where M~ Y x := {y : My = x} and vol n (M _1 a; n R") denotes the volume of 
M _1 xnIR™ in R™. In the following, we review some basic properties of multivariate 
truncated powers. For more detailed information about this function, the reader is 
referred to [8j [15] . 

For y = (j/x, . . . ,y s ) € R s and a function / defined on R s , we denote by D y f 
the directional derivative of / in the direction y, i.e., D y — ^2 S j = iyjDj, where 
Dj denotes a partial derivative with respect to the jth coordinate. The following 
differential formula was given in [8]. For y g M, we have D y T(-\M) = T(-\M\y). 
More generally, 

(2.3) D Y T{-\M) = T{-\M\Y) for Y C M 

where Dy := \[ y&r D v . 

Let W (M) denote the set consisting of those Y C M such that M\Y does not 
span R s and let D(M) denote the linear space of those infinitely differentiable 
complex- valued functions / on R s that satisfy the following system of linear partial 
differential equations: 

D Y f = 0, V re W (M). 
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Based on the definition of D(M), we can see that D(M') C D(M) when M' C M. 
It was also proved in [51 [13] that D(M) C P n _ s . 

Let c(M) be the union of all the sets cone(M \ Y), as Y runs over W{M), A 
connected component of cone°(M) \c{M) is called a, fundamental M-cone according 
to [12]. Then T(-\M) agrees with some homogeneous polynomial of degree n — s in 
D(M) on each fundamental M-cone. In fact, T(-\M) is generally continuous and 
positive on cone°(M). 

We now turn to box splines. The box spline B(-\M) associated with M is the 
distribution given by the rule [131 114] 

(2.4) / B(x\M)(j)(x)dx = / <j>{Mu)du, <t> G 2>{W). 

JR 3 -'[0,1)" 

According to ([2.4|1 . the support of 2?(-|M) is [[M]], and hence, we have (see [15] . 
pp.33) 

{j G Z s : B(x - j\M) 0} = Z s n (ar - [[M))). 
By taking = exp(— in (|2.4p . we obtain the Fourier transform of B(-\M) as 

- LJ - iC rrin 

For more detailed information about box splines, the reader is referred to [15) . 

3. Discrete truncated powers 

Let M be an s x n matrix with integer columns mi, . . . , m„ and suppose that 
[M] does not contain the origin. From the definition of t(-\M) given in Section 1, 
we have 

(3.1) Y, <P(a)t(a\M) = ^ ip(M/3), for any ip G 2?(K j. 

A comparison between (|3.1|) and (|2.ip shows that (|3.1[) is a discrete version of (|2.1I) . 
This observation was the motivation to designate t(-\M) a discrete truncated power. 
Also the identity 

(3.2) T(x\M) = Y t(a\M)B(x - a\M) 
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established in [10 j shows that these three functions, T(-\M), B(-\M) and t(-\M), 
are closely related. 

We denote by S the linear space of all complex functions on Z s . Given y G 7L S , 
the backward difference operator V y is defined by the rule V y / := / — /(• — y), 
where / G S. More generally, for an integer matrix Y , we define 

Vy := n v «- 

y£Y 

In [12], the following difference formula was given: 

V„t(-|M) = t(-\M \ y), for any y G M. 

More generally, 

V Y t(-\M) = t(-|M\y), for any Y C M. 
Given (9 = (0 1; . . . , S ) G (C \ 0) s , we set 

Me := {y G M : 6» y = 1}. 

Let 

A(M) := {6 G (C \ 0) s : span(M fi ) = IR S }. 
As pointed out in [11 , 9 G A(M) if and only if it has the form 

(3.3) 9 = cxp(27ria/detF), 

for some Y G 3§(M), and the vector a G Z s satisfying y T a = |detY|L, where 
L G Z s n [[y T )), i.e., is a lattice point in the parallelepiped determined by Y T . 
Here, 

m(M) = {Y C M : #y = s, span(y) = R s }. 

We let 

V(M) := {/eS:Vy/ = 0, for all Fefftf)}. 
In [11], Dahmen and Micchelli showed that a sequence / G V(M) if and only if 
it has the form /(a) = J2eeA(M) @ a Pe a e <^ S j where G D(Mg) for each 
G A(M). For a fundamental M-cone fi, we set 

u(n|Af) := z s n(n-[[M}}). 

Dahmen and Micchelli also proved the following result in [T2"] : 
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Theorem 3.1. ( [121 Thm.3.1]) Let M = {mi, . . . , m n } be an s x n integer matrix 
and suppose that M spans R s and the convex hull of M does not contain the origin. 
Then for any fundamental M-cone f2, there exists a unique element fn(-\M) G 
V(M) such that fn(-\M) agrees with t(-\M) on v(fl\M). Moreover, fn(-\M) has 
the following properties: for any x G Q such that 

v{x\M) n cone(M) = {0}, 

fn(-\M) is uniquely determined by 

fn(a\M) = 5 0a , a G v(x\M), 

and satisfies the relation 

n 

(3.4) f a (a\M) = {-l) n -'f a {-a-Y i m j \M), a e Z s . 

3 = 1 

We denote by 9^pg the complex sequence given by a H> 9 a pg(a), with a£Z s , 
where pe G P(R S ) and 9 G (C \ 0) s . Let E denote the linear space of all sequences 
of the form / = ^2 g 9^pg. It is easily seen that each / G E can be written uniquely 
in the form ^2 e 9^pg. Hence, we can define mappings Pq from E into P(R S ) as 
Pe ■■ f^Pe- 

Based on Theorem 13. 11 fn(ct\M) G V(M). So, it can be written in the form of 

(3.5) fn(a\M) = ]T 9 a Pe (a), a G Z s , 

eeA(M) 

where pg(a) G D(Mg). According to the definition of Pg, one has Pgfa,(a\M) = 
pg(a), for any 9 G A(M). Let e := (1, 1, . . . , 1) G Z s . As, for any y G M, e y = 1, we 
see that e G A(M). In fact, F e /o(-|M) is the polynomial part of f n (-\M). In [12] , 
Dahmen and Micchelli showed that the following equation holds: for x G SI, 

(3.6) T(x|M) = P e /o(x|M) + ^ D u P e /n( : r : |M)(-i) |u| J D"S(0|M)/ U !. 

0<\u\<.n — s 

From this equation, the following theorem was proved in |12j . 

Theorem 3.2. ( [T2l Prop. 5. 3]) Under the conditions of Theorem \3.1\ the leading 
homogeneous term of P e fn(-\M) agrees on f2 with T(-\M). 
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4. An explicit formula for discrete truncated powers 

The objective of this section is to present an explicit formula for discrete trun- 
cated powers t(-\M) which is the central result of the paper. Throughout the rest 
of the paper, we use f2 to denote any particular fundamental M-cone. According 
to Theorem 13. 1[ the discrete truncated power t(-\M) agrees, on v(Cl\M), with an 
element in V(M), which is denoted as fn(-\M). We use fl to denote the closure of 
ft. From Lemma 4.4 in [22], we have Q n 7L S C v(Q\M). Note that t(a\M) = 0, 



when a € I s \ cone(M). Hence, to present an explicit formula for t(-\M), we only 
need an explicit formula for /q(-|M). 

In the following theorem, we give the polynomial part of fn(-\M), i.e., P e /o(-|M). 

Theorem 4.1. Suppose that p^^n is the homogeneous polynomial part of degree 
n — s — fi of P e fo,(x\M). Under the conditions of Theorem \3.1[ when x £ Q, we 
have 



(4.1)p Mi n(x) 



'T(x\M), M = 0, 

-££0* (E\u\= k -3D"pMx)(~i) H D u B(0\M)/ul) , 

1 < [i < n — s. 



Proof. Rewriting of Equation (|3.6|) yields 

T(x\M) - D u P e f n {x\M)(-i)^D u B{0\M)/u\ 

(4.2) 0<|«|<n-s 

= P e fa(x\M), xen. 
According to (|4.2|) . po.n(x) agrees with T(x\M) on ft. Since both sides of (|4.2[l arc 
polynomials, we can rewrite (|4.2I) as 
(4.3) 

n—s 7i— s 

Po,a(x)- Yl D u J2p»M x )H) ]u1dU B(0\M)/u\ = J2p»Ax)> * G RS - 

0<|u|<n— s A*=0 /i=0 

A comparison between both sides of the equation above shows that 

PnM x ) = D u p jtQ (x)(-i) M D u B(0\M)/u\ \ , l<yt<n-s. 

3=0 \\u\=v,-j J 

□ 
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Based on p.5[) . to write down an explicit formula for t(-\M) on fi, we only 
need a formula for R&fn(-\M) for each fixed $ G A(M). To describe the formula 
conveniently, we suppose M\M$ = (mi, . . . , m K ), where k = #M— #M^. We select 
r as the least integer in the set {r : {d m ) r = 1 for all m G M \ M$}. (According 
to ()3.3|). the set is non-empty, since there is at least an integer, Yiye^M) 

|det(F)|, 

in it.) We use q® r (x), /i G Z + , to denote homogeneous polynomials satisfying the 
following condition: 

(4-4) <&.(»)= X! • • • D tT{x\Mo) ^ + ■ x £ Q. 

ii+-+u=n »=i ^ 

Throughout the rest of the paper, we set 

s 3 (x) := {-iy{x + Vx 2 + --- + {r-iy X r - 1 ) 

and M r := {mi,..., fh„}, where if m, G M#, then Sj = m*, otherwise, S, = rm,. 

Theorem 4.2. Suppose that is the homogeneous polynomial part of degree 
n — s — k — /i of P&fn(-\M). Under the conditions of Theorem \3.1[ for each fixed 
■d G A(M), we have 

At-l / 

<n(0 = C(-)-E E ^,n(0(-O |u| ^S(O|M r )/ U ! 
i=° \I«I=/*-j 

l</i<n — s — «;, 

where r is the least positive integer such that {§ m ) r = 1 holds for all m G M \ 

Combining Theorem 14.21 and qo ir (x) = T(x\ M#) Y[ w eM\M^ > x E CI, we 

can easily extend Theorem 13.21 as follows: 



Theorem 4.3. (|31j) Under the conditions of Theorem lff.il the leading part of 
Ptffn('\M) agrees with 

n-m n 

on fl. 



An alternative method for proving Theorem 14.31 is also given in [31) . 
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Remark 4.4. Since P^fsi(-\M) is a polynomial and S]nZ s C v(Q\M) (see Lemma 
4.4 in 22 J, Theorem O and Theorem S3] also hold on Q. 



Remark 4.5. In [5], Brion and Vergne give formulas for the lattice point enumerator 
of a convex rational polytope in terms of certain Todd differential operators. These 
are interesting connections to topology. The salient difference in approaches lies in 
our explicit use of splines. Although the formula presented in [5] might eventually 
be equivalent our, the analysis of the equivalence between our different formulas 
would be the content of another paper. 



Example 4.6. We consider the number of non-negative integer solutions of the 
linear equations 



3 2 10 
12 2 



M 

\X 4 / 



Set M = \ 2 2^ an< ^ k = (je^J ' ^ si m P ie computation shows that A(M) = 
{(1, 1)} U A 3 (M) U A 2 (M), where 



A 3 (M) = {(1, -1), (exp (2?ri/3), exp (27ri/3)), (exp (47ri/3), exp (47ri/3))}, 
A 2 (M) = {(-1,1), (exp (27ri/3), 1), (exp (2vri/3), -1), (exp (4ttz/3), 1), 

(exp (47ri/3), — 1), (i, —1), (exp (37ri/2), — 1), (exp (27ri/3), exp (57ri/3)), 

(exp (4ni/3), exp (ni/3))}. 



For 9 e A 3 (M), we have #Af e = 3 while #M e = 2 provided that 9 E A 2 (M). 
For the matrix M, there are 3 fundamental M-cones, i.e., Qi := cone°(Mi2), ^2 := 



cone°(M 23 ) and il 3 := cone (M 34 ), where M 12 = (jj J), M 23 = (J J j and M; 

1 

2 2 



34 
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Using Theorem 14.21 we obtain the explicit formula for t(-\M) as follows: when 
k€v{n u M), 



•exp((27rifci + 27rifc 2 )/3) + (-l) fe2 (fc 2 /24 + 5/48) + exp((47rifci + 4nik 2 )/3) 



+ exp(47riA: 1 /3)/18 + V3/l8i(-l) k2 exp(47ri/c 1 /3) + exp((2 7 rifc 1 + 57rifc 2 )/3)(l/24 
+V3/72i) + exp((47rifci + 7rifc 2 )/3)(l/24 - Vs/72i); 

when k e v(n 2 ,M), 

t(k\M) = 

fcifc 2 /18 - {k\ + fc 2 )/72 + 7fc 2 /72 + fei/18 + 23/108 + (-l) fe2 (fc 2 /24 + 5/48) 
+ exp((27rifci + 2mk 2 )/Z) ((3 + V3i)/108 • {2k x - k 2 ) + 19/216 + ^3/24^ 



+ exp((47rifci + 47rifc 2 )/3) ((3 - \/3i)/108 • (2fci - k 2 ) + 19/216 - Vs/24ij 
+(-l) fcl /16 + l/18exp(27rifci/3) - V3/18i exp(27rifci/3)(-l) fe2 + exp(47ri£;i/3)/18 
+V3/18i(-l) fc2 exp(47rifci/3) + 1/8 cxp(7rifc 1 /2)(-l) fe2 + l/8(-l) fc2 exp(37rifci/2) 
+(1/24 + Vsi/72) cxp((27riA:i + 57rifc 2 )/3) + (1/24 - ^^72) exp((4™fci + 7rifc 2 )/3); 

when k g w(ri 3 , M), 

t(k\M) = 

k\/2A + fci/4 + 47/144 + (-if 2 (fcj/12 + 1/4) + (-l) fcl /16 + exp(27rifci/3)/18 
-(V3/18)iexp(27rifci/3)(-l) fe2 + (1/18) exp(47rifci/3) + {V3/18)i exp(47rifci/3)(-l) fc2 
+(l/8)exp(7rifci/2)(-l)' £2 + (l/8)exp(37ri/ci/2)(-l) fc2 . 



In this section, we shall prove Theorem 14.21 The proof begins with several 
lemmas, designed to make the proof more readable. The main idea for proving 
Theorem l4.2l is to generalize Equation (|3.6[) . with the equation playing an important 



t{k\M) 





5. Proof of Theorem 14.21 
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role in obtaining the polynomial part for /n(-|M). This is the plan. Our discussion 
will be broken into three steps. 
First, we introduce a set. Let 

c (n,H) ■.= p| p| (n + b h h), 

heHb h e[os] 

where H is a finite set of real s- vectors. For example, if we set H := {1} and 
n := (0, +oo), then c(tt, H) = (1, +oo). We have 

Lemma 5.1. Suppose that H is a finite set of real s-vectors. Then for the funda- 
mental M-cone Q, the set c(f2, H) C £1 and the volume of c(17, H) in K. s is infinity. 

Proof. According to the definition of c(f2, H), one has 

c (n,H)=f] p| (n + b h h)c p| (n + b h h) = n. 

heH b h e[0,l] heH,b h =0 

Based on the definition of the fundamental M-cone, there exist real s-vectors 
01, ... ,5a,, oj e N, such that span{gi, . . . , g u } — R s and £1 — {a\g\ + ■ • • + a u g u : 
di > 0}. Hence, for any heH, there exist ai(h), . . . , a u (h) € R such that 
h = ai (h)gi -I 1- a w (h)g u . Then 

c(n,zr)=p| p U + 

fteffii),e[o ) i]ai,...,o u >o j=i 
Let a(/i) := max{|ai(/i)|, . . . , |a u (/i)|} and let a(H) := max{a(/i) : h e i/}. A 
simple observation is that dj + bhdi(h) > for any &/, € [0, 1], /i € if, provided that 
a, > a{H). Hence, we have 

S := {x : x = aih\ + • • • + a u h u , <ij € E, a, > a(-ff)} c c(0, i?). 

Note that the volume of S in R s is infinite. So, the volume of c(fi, iJ) in M s is also 
infinite. □ 

Second, we shall generalize Equation (|3.6[) . since it is of importance in obtaining 
the explicit formula for P e f$i(x\M). We recall some notations and definitions. Let 
§ be an arbitrary but a fixed vector in A(M). Recall that r is the least positive 
integer such that (d m ) r = 1 holds for all m £ M \ Ms- Let M r — {mi, . . . , fh n }, 
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where if mi G M$, then m« = m^, otherwise, fhi = rrm. Before extending Equation 
(|3.6|) . we first generalize (|3.2[) as follows: 

Proposition 5.2. 

(5.1) $~ a t{a\M)B(x-a\M r ) = ^ 0~ S"=i m ^T(x-J] m^M,.). 

a£Z s 0<ri ,...,r K <r z— 1 

Proof. Using the fact that i(a|M) = #{/? G : Mf3 = a}, we see that 
^ ?T Q i(a|M)B(a; - a|M r ) = ^ d~ Mfi B{x - M/3|M r ) = 

n 

(5.2) £ ^-^-^B^-^miAlMr). 
/3i,...,/3„ez+ i=l 

We can write /3j as 7;r + r^, where < r; < r, Jt,ri G Z + for i < k. Substituting 
Pi = JiT + Tii where i < K, into (|5.2I) and noting that $ rmi = l ) we obtain that 

E E E ^^-^S^-^m^ 

0<ri,...,r K <r- 7i,...,7 K eZ + / 9 K _ ( _ l5 ...,/3 n G^+ i=l 
k n 

- r 2J ^i7i - E TOi ^ 1 1 Mr ^ 

1—1 i — K-\-l 

K 

E $~^= im ^ ■ ^ B{x-M r p-^m iri \M r ) 

0<ri,...,r K <r /3eZ™ i=l 

K 

0<ri , . . .,r K <r i—1 

The last equation follows from X^ez" — M r /3|M r ) = T(x\M r ), which can 
be obtained directly from the definitions of B(-\M) and T(-\M). 
Now we arrive at 

K 

tf~ a t(a\M)B(x - a\M r ) = ^ 0-£"=i m * ri T(x - ^ m^lM;). 

ck£Z s 0<ri , . . .,r re <r i—1 

□ 



In fact, if we set § = e, then (|5.ip is reduced to (|3.2[) . So, as said above, (|5.ip can 
be considered as a generalization of (|3.2p . To simplify the term -d~ a t(a\M)B(x — 
a\M r ) in (l5~TT) . we need to study where 6 G A(M). 

Proposition 5.3. For any 6 G A(M), we Ziawe 0/i9 G A(M r ) and D(M e ) C 

D((M r ) e/ t). 
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Proof. The definition of M r shows that i) m = i?~ m = 1 for any m G M r . Hence, 
(M r )g = (M r )e/d for any e A(M). Also, rM e C (M r )e since m = 1 implies 
grm _ we k ave S p an ((M r ) e /^) = span((M r ) e ) = span(M fl ) = W which 

implies 9/d € A(M r ). To prove that D(M e ) C D((M r ) e /^), we only need to 
show that D(Af e ) c D((M r )g) since (M r ) e = (M r ) 8 /#. Select an / € D(Mg) and 
consider -D y /, where F = {y\, . . . , y w ) e ^ ((M r )e)- We claim that -Dy/ = 0. We 
set Y = (yi, . . . , y w ), where y rj = -jjj if jjj e M$ otherwise yj — y~j/r. We can see 
that Y € ^(Mg) and hence -Dy/ = c^Dyf = 0, where Co is a non-zero constant. 
So, we have / € D((M r )g), which implies that D(M e ) C D((M r )g). □ 

We now generalize Equation (|3.6p . Let := {(r — l)mj : 1 < i < n} and 
Q#,r(ar):= ^ $~ m,ri T(a; — m,r t |M r ). 
Then we have the following result. 

Lemma 5.4. When x E c(£l,H$), 

n — K — s 

(5.3) Q^, r (a;)=P^/n(x|M)+ £ C tl P^/ n (x|A/)(- i )l u l^ 11 S(.|M r )(0)/ U !. 

|«|=i 

Proof. By Proposition 15 . 21 we have 

Qtf.rfz) = ^ -i?" a t(Q!|M)B(a; - a|M r ). 
So, to this end, we only need to prove that 

n — K — s 

$- a t(a\M)B{x-a\M r ) = P#f n {x\M)+ £ ^"P fl /o(x|M)(-0 |u| ^ M 5(-|M r )(0)/ U !. 

aGZ s |u| = l 

The definition of c(H,,H^) shows that v(x\M r ) C t>(£l|M) provided that a; S 
c(0,JJ^). Noting that {j G Z s : B{x - j\M r ) ^ 0} = u(a;|M r ) (see J5], pp.33), we 
have 

(5.4) d~ a t{a\M)B{x - a\M r ) = d- a t{a\M)B(x - a\M r ), 

a£v(fl\M) aGZ 3 
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where x € c(0, Hence, 
(5.5) 

$- a t(a\M)B(x - a\M r ) = ^ $~ a t{a\M)B{x - a\M r ) 

aev(n\M) QGZ 3 

= 12 12 (0/V) a Pefn(a\M)B(x - a\M r ) + ]T P*/ n (a|M )B(aj - a\M r ) 
= ^ p $fn{oi\M)B{x - a|M r ) where x € c(Q,H#). 

The last equation follows from X) Q gz s p{ a )B(x — a|M r ) = for any p(a) € E(M r ) 
(see Proposition 5.2 in [12 ) and Ee G A(M)\tf (^/^"-Pe/nH^ ) G E(M r ), where 
E(M r ) is the space of functions of the form 

p(a) = 6a Pe( a )> Pe(a) e D((M r ) e ). 

8eA(M r )\e 

We now consider the sum 

(5.6) Y P-dMa\M)B{x - a\M r ), where x € c(Sl, H#). 

Let V(y) = P#f n (y\M)B(x - y\M r ). Then 
= / *(y)exp(-iyOdy 



P#fo(y\M)B(x - y\M r ) exp(-iy£)dy 

P&fn(x - t\M)B(t\M r ) exp(-i(x - t)£)dt 

= exp(-ix^) / P#fn{x-t\M)B(t\M r )exp(it£)dt 
= eM-ixt)P#h{-iD + x\M)B(-i\M r ). 

Taking into account that (see [5]) 

q(D)B(2na\M r ) = 0, where a £ I s \ and qE D(M r ), 

Poisson's summation formula converts the sum (|5.6[) into 

(5.7) Ptfn(-iD + x\M)B(-\M r )(0), 



1(5 
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since P$fn(- + x\M) G D(M$) for each fixed x and D(M$) C D(M r ). Expanding 
P#fn in a Taylor series for each fixed x and noting that _B(0|A/ r ) = 1, one has 

PM-iD + x\M)B(-\M r )(0) 
(5.8) = P*f a (x\M) + Yl D u Pvfn(x\M)(~i)^D u B(-\M r )(0)/u{, 

0< \u\<n — K — s 

where x G c(tt,H$). So, combining (j5"7T1) . (|5T3)| . (|53)l . (1531) . (pTTf and ([ST5| . we 

have 

n — k — s 

! 



Q*, r (aO = P*/n(a:|M) + £ ^"P^/o(a:|M)(-i)l u I^B(-|A/ J ,)(0)/ U 

|«|=i 



Finally, we present an explicit formula for Q^ >r [x). 
Lemma 5.5. When x G c(ft,H^), 



□ 



O^Ji H hire s — k. 



^ Cft + l)! 



Proof. When x £ c(fi, we can see that Q$ jr {x) is a polynomial of degree less 
than n — s + 1. Using the Taylor expansion, we obtain 

T(x - J2 nrriilMr) = T(x\M r ) + £ nD mi yT(x\M r ), 

2=1 j = l 3- i—1 

r-1 

where x G c(Q, Noting that ^ = o, we have 

Qa,r{x) 

K 

= 0-^ nmi T(x-Y,miri\M r ) 

0<ri ,...,tv <f i—1 



(5.9) 



E <r E <- nm< (r(x|M,,) + E^ f-E^) r(s|M r ) ) 

0<ri ) ...,r /s <r V j=l •'' V i=l / J 

0<ri,...,r„<r 0<jiH hj°„<n— s \i=l ) 

0<jiH \-j K <n-sO<Ti,...,r K <r \i=l / 
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Moreover, if there exists an index jh — 0, where 1 < h < K, then 
(5.10) 

k k / r— 1 \ 

E v-^=^ m > Ui-nDm^ = n E ^ rimi {-nD mi r = o, 

0<n,...,r re <r i=l i—1 \ri= / 

since E^=o ^"'""(-^-CmJ^ = E^=o ^ rh " lh = °- 

Recall that sj(x) = (-l) j (x + 2^x 2 H h (r - l^V"- 1 ). Based on (jO]) . (fSTTUjt 

and (12.31). we have 



o*,r(«)= E E (n ( " r< f. ]W) 

0<ri,...,r K <r 0<jH hj K <n-s \i=l / 

= e e fn^-^|p^V ( ^iM r ) 

~<ii+ - +j«<>«-s 0<n, ...,r*<r \i=l / 

= E ^ E - E (u^^^nxm 

k<31 + -~+3k<™-" 0<ri<r 0<r K <r \i=l Jl ' / 

0<jiH hi K <n-s-« i=l \0<ri<r KJt >' ) 

e ^.^rwn ^ , 

0<ji+— +i K <n-s-« i=l U * 

provided that x € c(f2, Htf). □ 

We now have all the ingredients for the proof of our main theorem. 

Proof of Theorem \4--2\ Since P§fn is a polynomial, it can be determined by its 
values on c(0, according to Lemma T5. II We recall (|4.4[) that, for /i € Z+, the 
homogeneous polynomials r (x) of degree n — s — k — ^ satisfy 

= E k-w^n^^.xeo. 

ii+-+i«=p »=i Ui ; ' 

Based on Lemma ESI we have 



0</x<n — s — k 

Note that P#/q(x|M) can be written in the form X)u=o _K ^ ) n( x )- By Lemma (5^3 
one has 

n — s — K n—K—s 

(5.11) E 4A X ) = P#fo(x\M) + Y, D U P#fo(x\M)(~i)^D u B(0\M r )/u\. 

M=0 |u|=l 
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Comparing the homogeneous polynomials on both sides of (|5.11j) . we arrive at 
Pin( x ) = <r(x)-Y,l E Dyl Q (x)(-i)WD"B(0\M r )/uA,t*>l. 

3=0 \\u\=n-j I 



□ 



6. Discrete truncated powers associated with special matrices 

In this section, we shall present more detailed information concerning t(-\M) for 
some particular matrices M. We first introduce some definitions. Let 

■Y k (M) := {Y C M : #Y = s + k, span(Y) = R s }. 

In this notation, &(M) = y (M). If gcd{|det(X)| : X € 3(Y)} = 1 for any 
Y G S^k(M), then M is called a k-prime matrix. In particular, when M is a 1- 
prime matrix, it is also called a pairwise relatively prime matrix. When s = 1, 
the fc-prime matrix has the property that any k integers in {mi, ... , m n } have no 
common factor. 

An explicit formula for t{-\M) is presented in [2] when s — 1 and M is a 1-prime 
matrix. 

Theorem 6.1. ([5]) Suppose that M — (ai, . . . , o„) € Z™ uraift ai, . . . , a„ pairwise 
relatively prime. Then 



n n „_ a 

t(a\M) = R^( ai ,...,a n ) + (-ir^2- E II 7' 

j=i ^ e a i=i^e fe=i,fc#j 

where R- a (ai,...,a n ) = -Res(F_ a (z), z = 1), and F^ a (z) = ^ 1 ^ z a 1 ). z .J_^ n ^ 1 _ z) ■ 
Based on Theorem 14.21 we can extend the theorem above to higher dimensions. 



Theorem 6.2. Under the conditions of Theorem \3.1[ 

h(a\M)=P e fn(a\M)+ E ^ | det L fl )l II ^^W^W, 

eeA(M)\e 1 ^ weM\M e 

provided that M is a 1-prime matrix, where P e fn(a\M) is presented in Theorem 

E3 
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Proof. Select a 9 £ A(M) \ e. To this end, we need to prove that #Mg = s. By 
(|3.3[) . 9 has the form 

9 = exp(27ria/| detY|), 
where the definitions of Y = (yi, . . . , y s ) and a are identical to those of (|3.3[) . Note 
that 9 V = 1 provided that y £ Y, i.e., Y C M g . We assert that Mg = Y. Otherwise, 
if there exists j/o € Af \ Y such that — 1, then we can select an s x s square 
matrix Y' C (yo,yi, . . . , y s ) such that yo € Y' and span(Y') = R s . Without loss of 
generality, we suppose Y' = {y\, . . . , y s _i, yo)- For each ?/ g Y', we have 9 V = 1. 
Hence, there is a /3 := (fit, . . . ,f3 s ) £ I s such that 9 is in the form of 

9 = cxp(27ri/3/|detY'|). 

Without loss of generality, we suppose that f3i ^ 0. Since < ol\ < \ det Y| — 1 
and < /3i < dct(Y') - 1, we have gcd(det(Y), det(Y')) > 1 due to a t /\ det Y| = 
/3i/| detY'|, which contradicts gcd(det(Y'), det(Y)) = 1. That is, Mg = Y. 

Note that Pgfa,(a\M) is a constant. This follows from the fact that the polyno- 
mial P e fn(a\M) £ D(M g ) C P . By TheoremgSl 

Pefa{a\M) = T(a\M e ) ]J j^^, 

weM\M(, 

when a £ Q. Then, we have 

fn(a\M) = P e fa{a\M) + ^ 9 a [] —L— T (a\M e ), a£QnZ s . 

9eA(M)\e w£M\M e 

It implies that 

h(a\M)=PJ n (a\M)+ ^ ]J _l_l conc(Me) (f7). 

0£A(M)\e 1 V n weM\M e 

□ 

We next show one example, as the application of Theorem 16.21 

Example 6.3. Consider the following partition function 

P{a,b}{n) = #{(%,y) £^ : ax + by = n,a,b £ Z + }, 

where a and b are relatively prime. We shall prove that P{ a ,b}( n ) = ^ — { b a ~ } — 
{ a b " } + 1, which is "the beautiful formula due to Popoviciu" [3], where {-} is the 
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fractional part of K b~ x and a -1 denote two integers satisfying b~ 1 b = 1 mod a, 
and o _1 a = 1 mod b. 

Whenn > 0,p {ajb} (n) = t(n\(a,b)). Note that T(x\(a, b)) = ^± and £> 1 _§(0|(a, &)) 
— | (a + b). By Theorem 16. 2\ when n > 0, P{ a ,h}(«) equals to 

n a + 6 1^ exp(^) 1^ exp(^) 
a& 2a6 a ^ 1 - cxp(-2z^fc) 6 i _ e xp(-^|^) ' 

k — J. " K — J. 

Based on the discrete Fourier analysis, one has (see [3], pp.144) 
1 j V 2a + a 2-l_exp(-^)' 

fc=l ^ v a I 

where t,a&Z. According to (JO]) , dHJ) can be reduced to ^ - { ^ } - { } + 1 . 

7. Exact volume of polytopes and the Ehrhart polynomial 

7.1. Exact volume of polytopes. A convex polytope P is the convex hull of 
a finite set of points in R™. In this section, we shall omit the qualifier "convex" 
since we confine our discussion to such polytopes. An integer polytope is a polytope 
whose vertices have integer coordinates. Similarly, a rational polytope is a polytope 
whose vertices have rational coordinates. The exact computation of the volume of 
P is an important and difficult problem which has close ties to various mathematical 
areas [TJ |31 [7J [25] • If P is a d-dimensional polytope in R" , then let vol„ (P) denote 
the d-dimensional volume of P in E™, i.e., the d-dimensional measure of P in M. n . 
Let Rp denote the affine space that is spanned by the vertex vectors of P. The 
lattice points in Rp form an Abelian group of rank d, i.e., Rp D 1 d is isomorphic 
to Z d . Hence, there exists an invertible affine linear transformation T : Rp — > R d 
satisfying T(Rp n Z") = 7L d . The relative volume of P, denoted as vol(P), is just 
the ci-dimensional volume of the image T(P) C R d . For more detailed information 
about the relative volume, the reader is referred to [55] . 

We next introduce a method for computing the relative volume of polytopes, 
which depends on the counting function for the integer points in a polytope. We 
consider the function of an integer-valued variable g that describes the number of 
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lattice points lying inside the dilated polytope gP: 

Lp(g) :=#(.gPf) Zn )- 

In [T5], Ehrhart inaugurated the study of the general properties of Lp(g). He 
proved that for an n-dimensional polytope P, Lp(g) is a polynomial in the positive 
integer variable g and that in fact 

(7.1) L P (g) = vol(P)g n + ±vol(dP)g n - 1 + ---+ X (P). 

Here, x(P) is the Euler characteristic of P and vol(dP) is the surface area of P 
normalized with respect to the sublattice on each face of P. Moreover, the leading 
coefficient of Lp(g) is the relative volume of P. Hence, if we obtain the leading 
coefficient of Lp(g), we can know the relative volume of P. In [I] and [7], the 
leading coefficient of Lp(g) was computed by the interpolation and the residue 
theorem respectively. We can even present a formula for the leading coefficient of 
Lp(g) using discrete truncated powers. Our result is 

Theorem 7.1. Suppose that PM,b '■= {x : A/x = b,x e M™} is an (n — s)- 

dimensional rational polytope. Here, M is an s x n integer matrix and b is an 
s-vector. Then vol(P M ,b) = Co • T(h\M), where C = gcd{| det(F)| : Y G @{M)}. 

Proof. To prove the theorem, we only need to prove that 

voi(p M , b ) - #{p-/ T )) n Z"} • T(b\M) 

since Cq = #{[[M T )) n Z' 1 } (see [30 ) . Since the dimension of Pjvf.b is n—s, we have 
b G cone°(M). By (|2.2p . T(h\M )/vol(PM,b) is a constant which is independent of 
b G cone°(M). For the integer matrix M, there exists an integer vector bo G 
cone°(M), such that Pjw.bo is an integer polytope. We now consider Lp Mbo (g). 
Through the definition of Lp M bQ (g) and t(-\M), we have 

Lp M , b(> (g) = t(gb \M), 
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with gho 6 cone°(M) when g G Z+. Based on Theorem 13.21 and Theorem 14. 3[ the 
leading term of Lp M bQ (g) is 

T(gb \M)l J2 8° b A = g n - s T(b \M) I £ 6*° 

\{9eA(M):M e =M} J \{9eA(M):M e =M} 

Since Lp MbQ {g) is a polynomial, we obtain that J2{eeA(M)-M g =M} ^ ffb ° ^ s a con " 
stant, i.e., 6 b " = 1. Hence, we have 

gisbo = ^ e A ( M ) . Me = M | = ^{[[m t )) n z n }. 

{eeA(A/):A/(,=A/} 

So, 

#{[[Af T )) n Z"} • T(bo|M) = vol(P M ,b ). 

It follows that 

T(b|M)/vol(P M , b ) = T(b |M)/vol(P M , bo ) = l/#{[[Af T )) fl Z"}, 
when b € cone°(M). Hence vol(P M , b ) = #{[[Af T )) n Z"} • T(b|M). □ 

Theorem 7.2. Suppose that P\ = {x : Ax < b,x G K" } is a n- dimensional 
polytope, where A is an sxn integer matrix and b is art s -vector. Then vol n (P b ) — 
T(b\M), where M = (A,E SXS ). 

Proof. Since P\ is an n-dimensional polytope in R™, the relative volume of P\ is 
equal to the volume of P\. Let P A:b = {(x,y) : Ax + E sxs y = b,x e M",y e 
Since E sxs C M, we have gcd{| det(Y)| : Y e ^(M)} = 1. If each component of 
b is rational, then P A ,h is a rational polytope. Note that vol n (F^) = vol(Pj) = 
vol(PA.b)- By TheoremOl vol n (P_£) = vol(P A ,h) = T(b\M) when P A>b is a ratio- 
nal polytope. Since both T(b\M) and vol„(P4 i b) are continuous at b and the real 
number can be approximated by rational numbers, one has vol„(P^) = T(b\M) for 
any b G cone" (M). □ 



How can T(-\M) be computed? In [24] . an efficient method for calculating the 
multivariate truncated power is presented. 
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Theorem 7.3. ([21]) Let M be an s x n matrix with columns mi, . . . , m n eZ s \0 
such that the origin is not contained in conv(M). For any Ai,...,A„ G R, and 

1 - 

(7.2) 2>|M) = _-^A 3 T(x|M\m J ). 

3=1 

Hence, one can compute T(x|M) according to the recurrence (I7.2|) . Combining 
Theorem 17.21 and Theorem 17.31 we can present an iterative method for computing 
vo\ n (P\). In fact, the computational complexity is 0{s n ~ s ). So, when ol := n— s is 
fixed, i.e., the dimension of the corresponding polytope is fixed, the run time of the 
algorithm based on recurrence (17.21) is a polynomial function of s. For a detailed 
description of the algorithm based on (|7.2[) , the reader is referred to the [19] and 

In the following, we shall present a simple proof for the volume formula for a 
polytope using this method, which is the central result in [26j . 
Let 

n„(.T) := {y e R" : Vi > 0, y 1 H h y { < x x H h x>, for all 1 < i < n} 

for arbitrary x := {x\, . . . , a; n ) with Xi > for all i. Let V^(a;) := vol(II„(a;)). The 
function V n (x) is a homogeneous polynomial of degree n in the variables Xi, . . . ,x n . 
This polynomial arises in a variety of mathematical fields, such as the calculation 
of probabilities derived from empirical distribution functions, the study of parking 
functions, and plane partitions. In [26], an explicit formula for V n {x) is presented 
using a probabilistic method. Based on (17.21) . we shall show a simple proof for the 
explicit formula for V n (x). 

Theorem 7.4. ([26 ) For each n = 1, 2, . . . , 



x , 



k£K„ kGK„ 

where K„ := {k G Z™ : J2l=i — 3 ^ 0T a ^ 1 < J < n — 1 and ^" =1 fe, = n}. 

Proof. Let b = (&i, . . . , b n ) T , where hi = X)/i=i a; ' i ' ^ denote the ith column 
of the n x n identity matrix, and let Li = X\ =i ifj. Let Mq = (L±, . . . , I n ) and 
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M = (Mo, E nxn ). Using these notations, we have 

n„(x) = {yeRl :M y<b}. 

By Theorem 17. 21 V n (x) — T(b\M). To this end, we present an explicit formula for 

T(b\M). Note that b = x x h + {x x + x 2 )h H h (ELi x h)I n - By Theorem[LS 

n 

T(b|M) = -{x x T{b\M \ h) + ■ ■ ■ + (V x h )T(b|M \ /„)). 

Note that 



b = x x h + (xi + x 2 )h H h (2J ^/i)^n 

i— 1 i n 

= cci/i H h x h )Ii-i + (2J ^h)-^ + afi+i-ft+i H 1- ( 2^ 

/i=l /i=l h=i+l 

i— 1 i n 

= xih H h Xh)h-i + (2J + 2^+1-^+1 H h ( ^ XfcKn- 

/i=l /i=l h=i+l 

For any fixed integer i, where 1 < i < n, we have T(h\M \ 7j U Ij) = since 
b ^ cone(Af \ Jj U ij). Hence, 

1 4-1 
T(b|M\Ji) - -(a: 1 T(b|M\J i UJ 1 ) + "- + (y)a! h )r(b|M\J i UJ i _i) 

+x i+1 (b|M \ Z, U + ( ^ x h )T(b\M \ h U /„)). 

h—i+l 

Moreover, for any i > j, we have 

j'-i j 
b = .Tiii H h (2J x h )Ij-\ + (^2 x h) I j + Xj+ilj+i H h 

h=l h=l 
i—1 i n 

( 51 affc)/j_l + ( ^ 3f/i)Ii + X i+ il l+ i H h( ^ X h)In- 

h=j + l h=j + l h=i+l 

Substitute b into T(b\M\(IiUlj)) and then expand T(b\M\ (IjU/,-)) by Theorem 
17.31 Continuing the process, we have 

V n {x) = T(b\M) = - ]T x ^ 
n\ z — ' 

where a n is the set of permutations of the set {1, 2, . . . , n}, x(i) — x(i\, . . . , i n ) — 

x(ii) ■ ■ -x(i n ) and x(i k ) = x 3k+x H Vx ik . Here, j k = max J k (when J k = 0, j k = 

0), and J k = Uj<fc-i ij<i k {h} f° r a nxe d permutation i = (ii, . . . ,«„). To this end, 
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we prove 

(7.3) £>(i) = e 

by induction. When n = 1, it is easy to prove that the conclusion holds. Suppose 
that the conclusion holds for n — m, i.e., 



x k 



A: 



We now consider the case where n = m + 1. Note that x m +i appears in x{ik) if and 
only if ik = m + 1. Hence, by (|7.4j) . we have 



E = E afea;fc ' 

where, if fc m+ i = 0, 

/ m \ / m \ /m + 1 

y/ci 1, . . . , fc m y V^i) • ■ • j 1/ \ ™ 

if fc m+ i = 1, then a k = (m + l)( kl ™ k J = ( kll .X\ m J = CT)- Noting 
km+i < 1, we have 



E «(«) = E 



Hence, (|T.3[) holds when n = m + 1. □ 

7.2. An explicit formula for Ehrhart polynomials. The explicit formula for 
Lp(g) is interesting. For any rational polytope, Ehrhart proved that Lp(g) is a 
quasipolynomial in g. Here a quasipolynomial is an expression of the form c n (g)t 9 + 
• • • + co(.g), where the Ci(g) are periodic functions in g. In (|7.ip . three coefficients of 
Lp(g) are presented. The other coefficients of Lp(g) have remained a mystery, even 
for a general lattice 3-simplex, until rather recently with the work of Pommersheim 
[25] in R 3 , Kantor and Khovanskii [23 in M 4 , Cappell and Shaneson 6: in W\ 
and Diaz and Robins [17] in R n . In the following theorem, the explicit formula for 
Lp(g) is presented in terms of multivariate truncated powers, where P is a rational 
polytope. 
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Theorem 7.5. Let {n — s)- dimensional polytopes P — {x £ K™ : Mx — b}, where 
M is an s x n integer matrix and b is an integer s -vector. We have 

n-s n-s-(#M-#M g ) 

I. run $>,n(b) 5 " 5 U Yl £ (6 h y P ln(b)g n - s -(* M -* M ^- j , 

j=0 9eA(M)\e j=0 

where £1 is a fundamental M-cone, pj_a(x) and p® n (x) are presented in Theorem 
\4-l\ and Theorem \4- e ^\ respectively. 

According to the definition of Lp(g) and t(-\M), we have Lp(g) = t(gh\M). By 
Theorem 14.11 Theorem 14.21 and the properties of multivariate truncate powers, the 
theorem follows. 
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